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The possibility of using the two-fold topological degeneracy of spin-1/2 chiral spin liquid states 
on the torus to construct quantum error correcting codes is investigated. It is shown that codes 
constructed using these states on finite periodic lattices do not meet the necessary and sufficient 
conditions for correcting even a single qubit error with perfect fidelity. However, for large enough 
lattice sizes these conditions are approximately satisfied, and the resulting codes may therefore be 
viewed as approximate quantum error correcting codes. 
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I. INTRODUCTION 



If one could be built, a quantum computer would be 
capable of solving certain computational problems much 
more efficiently than any classical computer, most no- 
tably factoring large integers into primes in polynomial 
time [p]| and searching unordered lists of N items in 
0(N 1 '^) queries In a quantum computer classical 
bits, which take the values or 1, are replaced by quan- 
tum bits, or qubits — two-level quantum systems whose 
Hilbert spaces are spanned by the orthonormal states |0) 
and |1). Unlike a classical bit, a qubit can therefore be 
placed in an arbitrary quantum superposition a|0)+/3|l). 
However, due to the coupling of this qubit with the out- 
side world, which may be small but which can never be 
reduced to zero, this state will eventually become entan- 
gled with its environment, losing its quantum coherence. 
Because maintaining this coherence is crucial for quan- 
tum computers to achieve their superiority over classical 
computers, the question of how to protect qubits from de- 
coherence has been central to the ongoing effort in quan- 
tum computing. 

One of the most surprising recent developments in 
quantum information theory has been the discovery of 
a scheme for fighting decoherence using what are called 
quantum error correcting codes P^|. A quantum error 
correcting code is a mapping from the Hilbert space of 
a single qubit, (or, possibly, more than one qubit), to a 
subspace of the Hilbert space of many physical qubits. 
The resulting many qubit state is then referred to as an 
encoded qubit. These encoded states are carefully de- 
signed so that if an error occurs, i.e., if a small number 
of the physical qubits become entangled with their en- 
vironment, certain measurements can be performed to 
determine which error has occurred and how it can be 
corrected without disturbing the quantum information 
stored in the encoded qubit. 

An important connection between quantum error cor- 
recting codes, many-body physics and topological quan- 
tum numbers was pointed out by Kitaev Q who con- 
structed a class of spin Hamiltonians realized on two- 
dimensional lattices with periodic boundary conditions. 



For the simplest of these models the physical qubits cor- 
respond to spin-1/2 particles located on lattice edges. Ki- 
taev showed that these models have degenerate ground 
states which are distinct, orthogonal states but which 
nevertheless cannot be distinguished by any local oper- 
ators. Instead, these degenerate states can only be dis- 
tinguished by global operators, e.g., operators which act 
on a set of spins which encircle a topologically nontrivial 
orbit of the torus formed by the periodic lattice. Kitaev 
showed that these 'topologically' degenerate states could 
be used as quantum error correcting codes, called toric 
codes, in which encoded qubits correspond to particular 
superpositions of these degenerate ground states. Be- 
cause these states can only be distinguished by topologi- 
cal quantum numbers, and not by any local observables, 
the environment has difficulty 'measuring' the encoded 
quantum information which is therefore, to some extent, 
protected from decoherence. 

This connection between quantum spin models with 
topologically degenerate ground states and quantum er- 
ror correcting codes provides motivation to revisit some 
of the topological degeneracies which can arise 'naturally' 
in certain condensed matter systems. In this paper, as 
a concrete example, the spin-1/2 chiral spin liquid states 
originally proposed by Kalmeyer and Laughlin [|| as pos- 
sible spin liquid ground states for frustrated spin-1/2 
quantum antifcrromagnets are considered in this context. 
It is known that these states possess a kind of topologi- 
cal order Q , not unlike the topological order of Kitaev's 
toric codes, which leads to nontrivial ground state degen- 
eracies on Riemann surfaces with genus 1 or greater. It 
should be emphasized that it is by no means clear that 
the results of this paper will be useful for constructing 
quantum error correction schemes for realistic quantum 
computers. Rather, the goal of the present work is to pro- 
vide some insight into the possible ways that the Hilbert 
space of an array of qubits can exhibit topological quan- 
tum numbers. 

The paper is organized as follows. In Section II the ba- 
sic physics of the chiral spin liquid states is reviewed. The 
case of finite Ni x N2 periodic lattices is considered and 
it is proven that, if properly constructed, the chiral spin 
liquid states realized on these lattices are exact singlet 
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states, generalizing a previous proof due to Laughlin that 
these states are singlets for TV x N periodic lattices with 
N even [0. In addition, it is shown by explicit construc- 
tion that these states possess a topological degeneracy 
on any periodic lattice, in agreement with [Q. In Section 
III the nature of this topological degeneracy is character- 
ized using the Lieb-Schultz-Mattis 'slow twist' operator 
and it is shown to be related to a topological decoupling 
of the Hilbert space of short-range valence-bond states 
on periodic lattices. In Section IV the general proper- 
ties of quantum error correcting codes are reviewed and 
compared with analogous properties of the topologically 
degenerate chiral spin liquid states, which are computed 
by variational Monte Carlo. It is shown that, unlike the 
toric codes, chiral spin liquid states do not satisfy the 
necessary and sufficient conditions to be exact quantum 
error correcting codes. However, as the lattice size in- 
creases, the violation of these conditions becomes weaker 
until, in the thermodynamic limit, they become satisfied. 
In this sense the topologically degenerate chiral spin liq- 
uid states may be viewed as approximate quantum error 
correcting codes on large enough lattices. Finally, Section 
V summarizes the results and conclusions of the paper. 



II. CHIRAL SPIN LIQUID WAVE FUNCTIONS 
ON FINITE PERIODIC LATTICES 

Consider the following spin-1/2 Hamiltonian realized 
on a two-dimensional square lattice, 

h = Jt J2 Si • Sj + j 2 s « ■ s j> w 

(id) {(id)) 

where J\, J 2 > and and denote nearest- 

neighbor and next-nearest-neighbor pairs of lattice sites, 
respectively. The lattice size is taken to be Ni x N2 
with Ni even and lattice spacing b. Periodic boundary 
conditions will be assumed throughout the paper. For 
concreteness (and future reference) the lattice is taken to 
lie in the xy plane, with lattice sites r = (tt-ix + n 2 y)b 
where n± and ni are integers. 

If Ji > and J2 — Hamiltonian (|I]) describes an un- 
frustrated two-dimensional spin-1/2 Heisenberg antifer- 
romagnet for which the ground state is known to possess 
long-range Neel order in the thermodynamic limit. In the 
opposite extreme J2 ^> J\ the two sublattices decouple, 
and each develops Neel order independently. It is gen- 
erally believed that over an intermediate range of J2 / J\ 
values the ground state is in a 'spin-peierls' phase with 
a locally observable broken translational symmetry, but 
there is no evidence that Hamiltonian (|l]) ever has a spin 
liquid ground state, i.e., a ground state with neither long- 
range Neel order nor any other locally observable broken 
translational symmetry, other than at zero temperature 
critical points. Nevertheless, in what follows Hamilto- 
nian (pi) will be used to introduce the chiral spin liquid 



states with the understanding that while these states al- 
most certainly do not describe the ground state of (|lj) in 
any parameter range, they may be eigenstates of an, as 
yet unknown, frustrated spin Hamiltonian. Fortunately, 
the topic of this paper - the relationship between chiral 
spin liquid states and quantum error correcting codes - 
involves properties of Hilbert space and does not depend 
on the Hamiltonian. 

Hamiltonian ([!]) can be viewed as describing a sys- 
tem of N interacting hard core bosons hopping on a 
square lattice where N = N1N2/2. In this description 
the bosons correspond to up spins moving in a down 
spin background with matrix elements J\ and J2 for 
nearest-neighbor and next-nearest-neighbor hopping, re- 
spectively. If the totally symmetric wave function de- 
scribing these bosons is <I>({ri}) then the corresponding 
spin state is 

|$)= Y, *({rm + A + 2 -'-S+ IU---D- (2) 

{ ri T 2 -r N } 

Because Ji and J 2 are both positive the effective hop- 
ping for the bosons is frustrated. As pointed out by 
Kalmeyer and Laughlin || this frustration can be viewed 
as being due to the presence of a magnetic field. To see 
this imagine that each bosons has charge q and moves 
in the presence of a magnetic field perpendicular to the 
plane of the lattice with field strength B — 4>o/b 2 where 
4>o = hc/q is the flux quantum. The corresponding vec- 
tor potential in the the Landau gauge is then A = —Byk 
and it can readily be shown that 

%2-k f Tj 

Ji = Ji exp — / A (21, (3) 

00 Jr, 

J 2 = - J 2 exp — / A • dl, (4) 

where r.j and Vj denote the starting and ending sites of 
the relevant hopping process and the line integrals are 
taken along straight lines connecting these sites. Thus 
a positive (frustrated) J2 corresponds to a negative (un- 
frustrated) J2 in the presence of a fictitious magnetic field 
of suitable strength. 

The sign of Ji can be changed without affecting J2 
by dividing the square lattice into A and B sublattices 
and rotating the spins on the A sublattice by 2ir radians 
about any fixed axis in spin space while leaving the B 
sublattice untouched. Under this sublattice rotation the 
boson wave function is transformed according to 

*({r,})= (^[e^+y^^}). (5) 

It is important for what follows to note that the spin 
wave function $ must satisfy periodic boundary condi- 
tions on an Ni x N 2 lattice in the x and y directions 
where Ni is even. Therefore for even values of N2 the 
transformed wave function \1/ must also satisfy periodic 
boundary conditions in the x and y directions, while for 



2 



odd values of N%, due to the sublattice mismatch, iff must 
satisfy periodic boundary conditions in the x directions 
and antiperiodic boundary conditions in the y direction. 

This mapping from a frustrated spin model to hard 
core bosons hopping on a lattice in the presence of a 
magnetic field inspired Kalmeyer and Laughlin to pro- 
pose a trial wave function based on the related system 
of interacting bosons moving in free space || . Following 
their work, imagine that no lattice is present and that 
the bosons move on a torus of length L\ — N x b in the 
x direction and L 2 = N 2 b in the y direction. Since the 
magnetic field is B = cf>o/b 2 = 2(hc/q)n, where n is the 
number density of bosons, the effective Landau level fill- 
ing fraction for the bosons is v = 1/2. In what follows 
the effective magnetic length Iq = (Hc/qB) 1 / 2 is taken as 
the natural length scale and set equal to 1, so that, for 
example, the lattice spacing is b = (2tt) 1 ^ 2 . 

A natural Ansatz for the ground state of this many- 
boson system with strong short-range repulsion is that 
the bosons condense into a v = 1/2 bosonic Laughlin 
state. The corresponding Laughlin wave function is com- 
pletely determined by the lowest Landau level constraint 
and the requirement that the wave function vanish as 
~ (zi — zj) 2 as two bosons approach one another. In the 
Landau gauge the Laughlin wave function for v = 1/2 
bosons on an L\ x L 2 torus can be written 

^({ Zi })l[e-ti/ 2 , (6) 

i 

where 

<&({*,}) = F(Z)HM<^ - (7) 

i<j 

Here $i(z|r) is the odd elliptic theta function ]l(J with 
t = 1L2/L1, Zi = Xi+iyi is the complex coordinate of the 
ith boson, Z = ^ z, is the center of mass coordinate, 
and 



F(Z) = e iKZ f[ tf a (n(Z -w v )/L x \t) 



(8) 



is the center of mass part of the wave function. 

The constants K, w% and u> 2 in (||) must be chosen so 
that iff satisfies the twisted boundary conditions 

*(rx + xii, r 2 , ■ • •) = e** 1 * (n,r 2 , • ■ ■), (9) 
¥(r x + yL 2 , r 2 , • • •) = e^e^^n, r 2 , • • •), (10) 

for each boson coordinate. Here 4>\ and 2 are two 
toroidal fluxes which characterize the x and y bound- 
ary conditions. The requirement that the wave function 
satisfy these boundary conditions leads to the following 
restrictions on K, w\ and W2 fl9[, 



e iKZ — e i( t>i 

e i2ir(w 1 +w 2 )/L 1 _ e i<j>2 e kL 2 



(11) 

(12) 



which can be satisfied in a variety of ways | 
K = 'coherent states' are used, for which 



F n {Z) = MAZ- W n )/L x \rf 



Here the 



(13) 



where 



W n 



-in 



Li, 



(14) 



for n = and 1. 

The two degenerate Laughlin states, iffo and ^i, cor- 
responding to n — and 1, are distinguished only by the 
difference in the center of mass parts of their wave func- 
tions, F n (Z). As shown by Haldanc Jll| this two-fold 
degeneracy is required for any translationally invariant 
system on a torus at v — 1/2. Note that although the 
states ^0 and ^1 span the two-dimensional Hilbert space 
of Laughlin states on the torus, they are not orthogonal. 
However, when L\ 3> £ 2 it can be shown that 



F n (Z) 




exp - 



TT{Z-W n -L 1 {m+\)f 
L1L2 



(15) 



This function is sharply peaked at the points Z = 
W n + L±/2 + mil for any integer m with peak widths 
~ \JL\L2- It follows that in the limit Li/L 2 — ► 00 the 
overlap between Fq and F\ vanishes and the states ^0 
and if>i become orthogonal. 

To use ^0 and ^1 as wave functions for hard core 
bosons on the Ni x N 2 periodic square lattice the bo- 
son coordinates are restricted to the lattice points r = 
(nix + n 2 y)6. For these lattice points the e~ tL2X fac- 
tor in ([h]) is identically 1 and the toroidal fluxes (j>\ and 
(j>2 then correspond to overall phases associated with the 
boundary conditions in the x and y directions. As shown 
above, for even and odd values of iV 2 the boson wave 
functions iff n are required to be, respectively, periodic 
and antiperiodic, in the y direction. Thus <f>2 = for 
even values of iV 2 and cf>2 = vr for odd values of 7V 2 , im- 
plying that the parameter W n appearing in F n (Z) is 



W n = 



nii/2 
(2n + l)£i/4 



A^ 2 even, 
N 2 odd. 



(16) 



Because N x is assumed to be even the wave functions 
iff n are always required to be periodic in the x direction 
implying that <pi — for all lattices. 

The spin states $„ obtained from iff n by undoing the 
sublattice rotation (|^) are referred to as chiral spin liquid 
states. These states break both time-reversal symmetry 
(T) and parity (P) under both of which $ -> $* [0. 
This broken symmetry, characterized by a nonvanishing 
chiral order parameter (a x ■ (er 2 x 03)), leads to a two- 
fold degeneracy which is clearly not topological because 
it can be identified by measuring a local order parame- 
ter, and therefore not potentially useful for constructing 
quantum error correcting codes. However, the additional 
degeneracy associated with the center of mass part of the 
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wave function is global in character and not associated 
with any local order parameter (at least in the thermo- 
dynamic limit, see Section IV). Using a field theoretic 
description, Wen has argued that the degeneracy of chi- 
ral spin liquids on a two-dimensional closed surface with 
genus g should be 2(2) 9 where the overall factor of 2 is 
due to the broken T and P symmetries and the factor 
of 2 9 is a measure the topological order j?J . The explicit 
construction of the states $o and $1 on periodic lattices 
is therefore consistent with Wen's result for g = 1. 

The requirement that the spin states $ n be singlet 
states is equivalent to the requirement S~ |$ rl ) = which 
in turn implies that 



^'*«(ri, 



P2, •••) = 0, 



(17) 



where the primed sum denotes a sum over lattice sites on 
the torus. If the toroidal fluxes are chosen so that $„ is 
periodic in both the x and y directions when evaluated 
on lattice sites then it can be shown that $ n satisfies dlT 
by using the singlet sum rule derived by Laughlin 



-M 2 /2 



0. 



(18) 



in z\. The chiral spin liquid states are therefore singlets 
for any N\ x N 2 lattice where N\ is even, provided the 
toroidal fluxes have been chosen, as they have been here, 
to ensure that the spin wave function is periodic in both 
the x and y directions. 

The nature of the topological degeneracy of the chiral 
spin liquid states depends on whether N 2 is even or odd. 
To understand this distinction consider the translation 
operator T x which translates each boson by one lattice 
vector in the x direction. Under this operator the rela- 
tive part of the wave function $„ is unaffected and only 
the center of mass coordinate is shifted according to 



Z + Nb- 



NiN 2 b 



N 2 L! 



(23) 



For odd values of N 2 this implies that the center of mass 
is shifted through a half-odd integer multiple of L\ and, 
due to the periodic boundary conditions, this is equiva- 
lent to a net shift of the center of mass by Li/2. Thus, 
in obvious notation, 



T x F n (Z) = F n (Z + L x /2), 
from which it follows that 



(24) 



where for lattice sites r = (nix + n 2 y)b, 
G(r) = (_iy' l i™2+«i+n2+i 



(19) 



and f(z) is any polynomial in z. Note that in order for 
( |l8| ) to be satisfied it is necessary to sum over all lattice 
points on the infinite two-dimensional plane. 

Following Laughlin j|] the sum rule ( |l8| ) can be applied 
to the chiral spin liquid wave functions for finite N% x N 2 
periodic lattices by first exploiting the periodicity of $ in 
the x and y directions to extend the summation in Jl7j ) 
to the entire lattice, 

W / \ , 2ttNiN 2 

\ $„(r l5 r 2 ,---) = hm _ \ $ n (n, r 2 , • • •), 



|ri|<fl 



e ib(x+ V )/2 e -y^/2 



-G(r)e 



^ 2 /4 e -M 2 /4 



T x $ n = $1 



and 



(25) 



This implies that for odd values of N2 the spin liquid 
states $0 and $1 break translational symmetry. In con- 
trast, for even values of N 2 , the center of mass coordi- 
nate is shifted through an integer multiple of L\ and the 
translation operator T x has no effect, 



T,$ = $0 



and 



(26) 



Finally, because N\ is even, for both even and odd values 
of N 2 



(20) 



T y <f> = $ and T y ^ = $ x , 



(27) 



and then using the following identity which holds for all 
lattice points, 



(21) 



where T y is the translation operator which translates each 
boson by one lattice vector in the y direction. 



III. CONNECTION TO LIEB-SCHULTZ-MATTIS 
OPERATOR AND VALENCE-BOND TOPOLOGY 



to show that 
53 $n(ri,r 2 ,---) 

|ri|<H 

V G^e^F^Z)^, z 2 , ■ ■ .)e-l^l 2 / 4 J] e- yV2 - 



\ri\<R 



(22) 



The topological degeneracy of the chiral spin liq- 
uid states can be elucidated further by introducing 
Affleck's two-dimensional generalization of the Lieb- 
Schultz-Mattis slow twist operator p^ , p^ |, 



Ulsm = cxp ( ij^^ 



(28) 



In the limit R — > 00 the summation on the right hand 
side of ( p2] ) vanishes due to the sum rule ( |l8| ) and the 
fact that the function e Zl / 4 F n (Z)?p(zi, z 2 , • • •) is analytic 



where the primed sum denotes a sum over lattice points 
on the torus. The usefulness of this operator derives 
partly from the fact that when Ni 3> N 2 , for any singlet 
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state | Sing) and any rotationally invariant spin Hamilto- 
nian such as ([!]) which only includes short-range interac- 
tions, it can be shown that |l3|]l4| ] 

(Sing| (u LSM HUi SM - H) |Sing) ~ O (j^j , (29) 

where J is a measure of the typical magnetic interaction 
strength. If, as is supposed to be the case here, H has de- 
generate singlet ground states separated by a gap from all 
excited states, then ( p9| ) implies that in the Nx/N 2 — * oo 
limit Ulsm maps states in the finite dimensional Hilbcrt 
space spanned by these states into one another. 

The Lieb-Schultz-Mattis slow twist operator can be re- 
cast in bosonic language as 



of the topological degeneracy of chiral spin liquids in 
terms of the valence-bond state basis. This basis con- 
sists of states in which spins are singlet correlated in 
pairs, which are said to be connected by valence bonds. 
While any singlet state may be represented as a linear 
superposition of valence-bond states, it is reasonable to 
assume that any singlet state, such as the chiral spin liq- 
uid states, in which the spin-spin correlation function de- 
cays exponentially with distance ||[| can by represented 
as a superposition of short-range valence-bond states. A 
short-range valence-bond state is a valence-bond state 
containing only bonds with lengths less than a specified 
length, or bonds with a distribution of lengths which falls 
off exponentially for long bonds. 



Ul 



sm = exp 




exp 



■ 27T Y 

i — X 
Li 



(30) 



where X is the x coordinate of the center of mass. Due to 
the periodic boundary conditions there is some freedom 
in labeling the lattice sites on the torus, and in order to 
precisely define U lsm it is necessary to choose a partic- 
ular labeling scheme. Here it will be assumed that the 
primed sum in ( |30| ) is over lattice sites r = (nix + n 2 y)b 
where n% = -N x /2 + 1, • • • , Ay 2 and n 2 = 1, • • • , AT 2 . 
For this choice l2' r x = ^2-^l/(2&) and 



Ulsm = (— i) 



I - 27r v 

exp [ i— X 



(31) 



As shown in Section II, when L\ 3> Li the center of 
mass part of the chiral spin liquid wave functions <E>„ be- 
comes sharply peaked for Z = W n + L\/2 + mL\ for any 
integer m. Therefore in this limit 



rxp[ij^X ) F.AZ) 



«'x P (i^W„ ) FjZ). {■■YD 



Combining ( |3l| ) and ( |32| ) and using (16) one finds that 
in the N1/N2 — * 00 limit the states $0 and $1 become 
eigenstates of Ulsm with eigenvalues ±1, 

{JVo+2 
(-1)"(-1) — N 2 even, 
(33) 
(-l)"(-l)^-$„ N 2 odd. 

Figure ^ shows the results of a numerical variational 
Monte Carlo computation of the real part of the ex- 
pectation values (& n \ULSM\$n) for n = and 1 and 
N 2 = 3 and 4 plotted vs 1/N\. The results clearly show 
that Re(<b n \ULSM\&n) — > ±1 in the Ai — > 00 limit for 
fixed N 2 . Because Ulsm is a unitary operator it fol- 
lows that $0 and $1 become eigenstates of Ulsm in the 
iV^A^ - ► 00 limit with eigenvalues ±1, consistent with 

The fact that $0 and <&i become eigenstates of Ulsm 
in the N\/N 2 — > 00 limit leads to an appealing picture 
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FIG. 1. Real part of the expectation value of the 
Lieb-Schultz-Mattis slow twist operator in the chiral spin liq- 
uid states $0 and $1 for Ni x 3 lattices (top) and Ni x 4 
lattices (bottom) plotted vs 1/N\. Statistical error bars are 
smaller than symbol sizes. 



The requirement that valence bonds must connect two 
sites, and only one bond may be attached to each site, 
gives rise to a topological decoupling of the space of short- 
range valence-bond states ]l^-18[]. Figure | shows four 
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short-range valence-bond states, two on a 6 x 3 lattice and 
two on a 6 x 4 lattice. In this figure solid lines connecting 
pairs of lattice sites represent valence bonds. In each of 
these states the x projection of the length of each bond 
does not exceed 2b and so it is possible to unambiguously 
determine the way in which a given bond 'wraps' around 
the periodic boundary condition in the x direction (it is 
in this sense that these states are short-range valence- 
bond states). For each of these states 6 vertical dashed 
lines are shown which 'slice' the gaps between each ver- 
tical line of lattice sites. The parity (o = odd, e = even) 
of the number of bonds cut by these dashed lines is then 
shown below each line. 
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FIG. 2. Four short-range valence-bond configurations illus- 
trating the topological quantum numbers responsible for the 
degeneracy of chiral spin liquid states on periodic lattices. 
Dashed lines are drawn through horizontal gaps in these con- 
figurations with the parity of the number of valence bonds 
crossed by each line shown below (e=even, o = odd). The up- 
per two configurations on 6 x 3 lattices indicate the generic be- 
havior for odd width lattices in which an alternating even-odd 
or odd-even pattern appears. The lower two configurations on 
6x4 lattices indicate the generic behavior for even width lat- 
tices in which the gap parities are either all odd or all even. 
The configurations on the left and right contribute, respec- 
tively, to the states $o and $i, in a sense described in the 
text. 

For N2 — 3, or any odd value of N2, an alternating 
even-odd pattern invariably appears plf . Short-range 
valence-bond states then fall into two distinct classes, 
which can be referred to as even-odd and odd-even, cor- 
responding to the two N2 = 3 configurations shown in 
Fig. ||. For N2 — 4, or any even value of N2, all gaps 
have the same parity and, again, there are two possibili- 
ties, either each gap has odd parity or each gap has even 
parity, corresponding to the two N2 — 4 states shown in 
Fig. ^|. It is convenient to define a topological quantum 
number, the gap parity, of a given short-range valence- 
bond state \a) to be (— l) 7a where "f a is the number of 
bonds which cross the gap between the line of lattice 



points with x — L\/2 and those with x = 1 — Li/2, i.e., 
those bonds which cross the discontinuity in x due to 
the periodic boundary conditions using the site labeling 
scheme introduced above. According to this definition, 
for the two configurations at the top of Fig. || j a = 2 and 
1, and the gap parities are +1 and —1, while for the two 
configurations at the bottom of Fig. || j a = 1 and 2, and 
the gap parities are —1 and +1. 

In |l8| it was shown that if N± ^> N2 for any short- 
range valence-bond state |a) 



U LSM \a) ~ (~1) 7 ». 



(34) 



The appearance of the gap parity, (~l) 7a , in (B4[) is due 
to the minus sign obtained whenever a spin- 1/2 particle 
is rotated through 2tt radians about any axis. If a given 
short-range valence-bond state is acted on by Ulsm then, 
if Ni 3> N2, for most valence bonds in |a) the two spins 
forming the bond are rotated by approximately the same 
amount. These valence bonds are therefore only weakly 
affected by the slow twist operator. However, for those 7 a 
bonds which cross the discontinuity in x due to the peri- 
odic boundary conditions the operator Ulsm rotates one 
spin by approximately 27r radians while the other spin is, 
again approximately, not rotated at all. Therefore, while 
these bonds also remain approximately in singlet states, 
they each contribute a factor of —1 to ([mJ) because only 
one spin has been rotated through 2-7T radians. 

According to (|3^), when N\ ^> N2 the states $0 and 
<&i become eigenstates of Ulsm with eigenvalues ±1. It 
is therefore plausible to assume that these states can 
be represented as linear superpositions of those short- 
range valence-bond states which also become eigenstates 
of Ulsm in this limit with the same eigenvalues, i.e., 
those states whose gap parities are equal to the corre- 
sponding eigenvalues given in (33). The gap parity can 
then be viewed as the topological quantum number which 
distinguishes between the states $0 an d $1 in this limit. 
Note that the alternating even-odd or odd-even patterns 
which appear in the gap parities for odd values of N 2 , 
and the uniform gap parities, either all even or all odd, 
which appear for even values of N%, are consistent with 
the symmetry properties of <3?o and $1 under the trans- 
lation operators T x derived in Sec. II. For more details 
on the connection between the Lieb-Schultz-Mattis slow 
twist operator and the topological decoupling of short- 
range valence-bond states see pl| . 



IV. CONNECTION TO QUANTUM ERROR 
CORRECTING CODES 



In the previous section it was shown that, in a sense 
which becomes precise in the limit N\ 3> N2, the topo- 
logical quantum number distinguishing degenerate chi- 
ral spin liquid states is the gap parity. This topological 
quantum number is similar to that of Kitaev's toric code 
in that it appears to be necessary to measure a global 



G 



property of the system, using, for example, the Lieb- 
Schultz-Mattis slow twist operator, in order to determine 
its value. Motivated by this similarity between chiral spin 
liquids and toric codes, it is natural to ask whether, or 
to what extent, the topologically degenerate chiral spin 
liquid states on finite lattices can be viewed as quantum 
error correcting codes. 

A quantum error correcting code for a single qubit is 
a mapping of the form, |0) — > |0l) and |1) — > where 
the states |0l) and are made up of several physical 
qubits. If the encoded qubit is placed in a pure state 
\Tl) = ci|0l) + then the initial density matrix de- 

scribing the state is po — |Tl)(Tl|. After the physical 
qubits making up the encoded qubit interact with their 
environment the most general effect on the density ma- 
trix is 



(?T.\<T?.\r, 



o, 



(41) 



Po 



J2E aPo El=p E = J2Ea\T L )(T L \El (35) 



with the constraint ^ a ElE a = 1 where the operators E a 
are referred to as error operators. In order to be able to 
return the encoded qubit to its original pure state there 
must exist a recovery operation which satisfies 



PE -> RaPsRl 



Po 



|Ti)(T z 



(36) 



again with the constraint Yla^a^-a = 1- The necessary 
and sufficient conditions for such a recovery operation to 
exist are MM 



(0 L \AlA b \0 L )={l L \AlA b \l L ), 
(0 L \A{A b \l L )=0, 



(37) 
(38) 



where the set of operators {A a } form a linear basis for 
the error operators, i.e., every error operator can be ex- 
panded as E a — ^2 b X ab A b . For example, for a code ca- 
pable of correcting only single qubit error one may take 
the basis {A a } to consist of the identity operator and all 
Pauli matrices acting on individual physical qubits. 

The chiral spin liquid states <&o and <&i are not orthog- 
onal, except in the limit N1/N2 — > 00. However, on any 
finite lattice it is possible to orthogonalize them and use 
the resulting states as a quantum code where the spin- 
1/2 particles located at lattice sites correspond to the 
physical qubits. For example, 



\0l) 



l*o), 



(*ol*i)l*o) 



(l-K$o|$i>| 2 ) 



1/2 



(39) 
(40) 



The question to be addressed is then, to what extent do 
these states satisfy the criteria ( |37| ) and (^|) for being 
quantum error correcting codes? 

Because $0 and <&i are singlets it is possible to sim- 
plify d37j ) and ( |38"| ) considerably for the case of single 
qubit errors by noting that an arbitrary encoded qubit 
|Tl) = a|0i) + /3| 1 x) is also a singlet, implying that 



for all lattice sites where a = x,y, or z, thus ensuring 
that ( |37| ) and d38| ) are satisfied for A a = 0% and A b = 1. 
Likewise, 

(r L \*«aP.\y L ) = s a0 {r L \^y T .\r L ), (42) 

for all lattice sites and r j . The conditions for a singlet 
state to be a quantum error correcting code capable of 
correcting a single qubit error can then be shown to be 
equivalent to the requirement that 



(T L \cj*ain L )=C r 



(43) 



for all states Tl, i.e., the spin-spin correlation functions 
must be identical for any encoded state. 

Because the spin-spin correlation function decays 
rapidly with distance in the chiral spin liquid states 
the largest violation of ( (43| ) is likely to occur for nearest- 
neighbor spin correlations. Consider these correlations 
for |T L ) = |$ ) and \T L ) = |$i). For odd values of N 2 
there is a broken translation symmetry in the x direction, 
T x $o = $1 and T x <&i = $01 and, for any lattice site ro, 

($n|« +6 al$n) = ($„ l< +26x< +6a+26i I *«> ( 44 ) 

and 

(*o|«+fcal*0> = (*lK+6x<+ba+6xl $ l)- ( 4 5) 

Here, and in what follows, a = x, y. For even values 
of N2 there is no broken translation symmetry in the x 
direction, T x $ n = $ n , and 

($»K<+6&I*») = < $ nK+fcx<+6a+6xl*n). (46) 

For both even and odd values of N2 there is no broken 
translation symmetry in the y direction, T y Q n = $„, and 
so in both cases 

< $ n|« +6a |*«) = (*nK+6y<+6a+6yl*n)- ( 4 ^) 

Finally, for odd values of N2 the fact that the chiral spin 
liquid states are symmetric under PT, the product of 
parity and time-reversal fll2|, implies that 



($oki a. 



Ol^ro^ro+fcyl^O 



H^ro^ro+by 



(48) 



Though it is not possible to compute these correlation 
functions analytically, it is straightforward to compute 
them numerically using standard variational Monte Carlo 
techniques, and the results of such calculations for vari- 
ous lattice sizes are given in Table |. In this Table the site 
ro is taken to be the origin (ro = Ox + Oy) and for each 
lattice size the correlation functions (^nl^ro^ro+ft*!^*) 
and (3?n|o'ro <J ro+6yl ( J ) » l ) are gi yen for n = and 1. Using 
the symmetry properties derived above, these correlation 
functions can be used to determine all nearest-neighbor 
spin-spin correlation functions for $ and <$>\. (Note that 
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for odd values of N2 there is some redundancy in the Ta- 
ble, since ($o|« n+by l$o) = <$i|« +6 yl^i)0 




FIG. 3. Patterns formed by nearest-neighbor spin-spin cor- 
relation functions (^nlo-^oy, |^n) for the topologically degen- 
erate chiral spin liquid states $o and $1 on 6 x 3 and 6x4 
lattices. Thicker lines correspond, qualitatively, to larger val- 
ues of — ($ n |o"r 4 o"rj | &n) (thicknesses of the lines are exagger- 
ated for clarity). For odd values of N2 the broken trans- 
lation symmetry is observable. For even values of N2 the 
nearest-neighbor spin-spin correlation functions are different 
in the two states. The ability to distinguish between $0 and 
$1 by measuring operators consisting of only two Pauli matri- 
ces indicates that although the underlying distinction between 
them is topological, as depicted in Fig. these states are not 
quantum error correcting codes on finite lattices. 

As can be seen in Table 0, on finite lattices the nearest- 
neighbor spin-spin correlation functions are not identical 
for $0 and $1, thus violating Therefore, on these 

lattices, the topologically degenerate chiral spin liquid 
states are not exact quantum error correcting codes, even 
for single qubit errors. While it is true that with increas- 
ing lattice size the difference between correlation func- 
tions in $0 and $1 becomes smaller, until it is no longer 
possible to distinguish between them due to the statisti- 
cal error bars of the Monte Carlo simulation, given the 
clear violation of ( |43| ) for lattices sizes as large as 8 x 6 it 
is unlikely that these correlation functions ever become 
exactly equal to one another on any finite lattice. Rather, 
it is more plausible that they approach each other expo- 
nentially as the system size, in particular iVjj, increases, 
though no proof of this has been found. 

The distinction between the states $0 and <E>i can be 
seen clearly in the patterns formed by the values of the 
nearest-neighbor spin-spin correlation functions. These 
patterns are shown for 6x3 and 6x4 lattices in Fig. [|. If 
the topologically degenerate chiral spin liquid states did 
provide exact quantum error correcting codes for single 
qubit errors then these patterns would be identical for a 
given lattice size. Figure together with Table | show 



clearly that despite the fact that the underlying distinc- 
tion between the states $0 and $1 is topological, as illus- 
trated in Fig. |^, on finite lattices the difference between 
them can still be measured locally using just two Pauli 
matrices. However, as stated above, the nearest-neighbor 
spin-spin correlation functions rapidly become effectively 
indistinguishable for these two states as the lattice size 
increases, as do, plausibly, all the correlation functions 
appearing in (f43|). In this sense the topologically degen- 
erate chiral spin liquid states on sufficiently large lattices 
may be viewed as approximate quantum error correcting 
codes. 



V. CONCLUSIONS 

In this paper the chiral spin liquid states first intro- 
duced by Kalmeyer and Laughlin as possible ground 
states for frustrated spin-1/2 antiferromagnets have been 
analyzed from the point of view of their connection to 
quantum error correcting codes. Explicit wave functions 
were constructed for the two topologically degenerate chi- 
ral spin liquid states on finite periodic N\ x N2 lattices 
with Ni even and it was proven that, if properly con- 
structed, these states are exact singlets for any such lat- 
tice. It was also shown that, in a sense which becomes 
precise when N\ 3> N2, the property characterizing the 
topological degeneracy is the gap parity — a topolog- 
ical quantum number associated with the short-range 
valence-bond state basis. However, despite the fact that, 
like Kitaev's toric codes, the degenerate chiral spin liquid 
states are distinguished by a topological quantum num- 
ber, these states are not perfectly indistinguishable when 
measured with local operators, except in the thermody- 
namic limit. Thus, on finite periodic lattices, these states 
do not satisfy the criteria (^) and (|38|), and so are not 
exact quantum error correcting codes — any error cor- 
rection scheme using chiral spin liquid states would not 
be able to recover even a single qubit error with perfect 
fidelity. However, the distinction between these states, 
as measured by local operators, rapidly becomes unob- 
servable as the lattice size increases. Therefore, on large 
enough lattices, the topologically degenerate chiral spin 
liquid states may be viewed as approximate quantum er- 
ror correcting codes. 
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TABLE I. 


Nearest-neighbor spin-spin 


correlation functions in 


the states <£>o and $i for different lattice sizes. 


Lattice Size 


<$o|« +M |$ ) 




X^lFroOro + Ml^l/ 




N x x N 2 










4x2 


-0.173(2) 


-0.946(2) 


-0.455(2) 


0.273(5) 


4x4 


-0.247(2) 


-0.246(3) 


-0.230(2) 


-0.376(3) 


4x6 


-0.216(2) 


-0.312(2) 


-0.217(2) 


-0.301(3) 


4x8 


-0.210(2) 


-0.306(3) 


-0.210(2) 


-0.307(3) 


6x2 


-0.176(2) 


-0.944(2) 


-0.467(2) 


0.322(5) 


6x4 


-0.311(2) 


-0.216(3) 


-0.279(2) 


-0.376(3) 


6x6 


-0.298(2) 


-0.300(3) 


-0.302(2) 


-0.281(3) 


6x8 


-0.303(2) 


-0.289(3) 


-0.302(2) 


-0.290(3) 


8x2 


-0.175(2) 


-0.944(2) 


-0.464(2) 


0.335(5) 


8x4 


-0.306(2) 


-0.210(3) 


-0.275(2) 


-0.382(3) 


8x6 


-0.290(2) 


-0.303(3) 


-0.292(2) 


-0.281(3) 


8x8 


-0.291(2) 


-0.291(2) 


-0.290(2) 


-0.293(3) 


4x3 


-0.230(2) 


-0.241(3) 


-0.301(2) 


-0.241(3) 


4x5 


-0.229(2) 


-0.301(3) 


-0.221(2) 


-0.301(3) 


4x7 


-0.213(2) 


-0.305(3) 


-0.213(2) 


-0.305(3) 


4x9 


-0.209(2) 


-0.306(3) 


-0.209(2) 


-0.306(3) 


6x3 


-0.334(2) 


-0.239(3) 


-0.257(2) 


-0.239(3) 


6x5 


-0.280(2) 


-0.290(3) 


-0.292(2) 


-0.290(3) 


6x7 


-0.283(2) 


-0.294(3) 


-0.281(2) 


-0.294(3) 


6x9 


-0.281(2) 


-0.293(3) 


-0.281(2) 


-0.293(3) 


8x3 


-0.258(2) 


-0.239(3) 


-0.336(2) 


-0.239(3) 


8x5 


-0.298(2) 


-0.290(3) 


-0.293(2) 


-0.290(3) 


8x7 


-0.291(2) 


-0.290(3) 


-0.292(2) 


-0.290(3) 


8x9 


-0.290(2) 


-0.291(3) 


-0.290(2) 


-0.291(3) 







